A Stochastic Model of Human Gait Dynamics 
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We present a stochastic model of gait rhythm dynamics, based on transitions between different 
"neural centers", that reproduces distinctive statistical properties of normal human walking. By 
tuning one model parameter, the hopping range, the model can describe alterations in gait dynamics 
from childhood to adulthood — including a decrease in the correlation and volatility exponents with 
maturation. The model also generates time series with multifractal spectra whose broadness depends 
only on this parameter. 



PACS numbers: 05.40.-a,87.15.Aa,87.90.+y 

Many physical and physiological processes exhibit 
complex fluctuations characterized by scaling laws , 
some with monofractal |2| and others with multifrac- 
tal behavior For example, normal gait dynamics of 
adults are monofractal B, while healthy heartbeat fluc- 
tuations are multifractal^ . The origin of these complex 
fluctuations and the factors that contribute to differences 
in their behavior are largely unknown. Here we develop 
a physiologically-motivated model that may be helpful 
in explaining some of the properties that contribute to 
complex dynamics. We focus on one class of signals — 
the time series of the inter-stride-intervals (ISI) between 
successive strides in human gait. 

Walking is a voluntary process, but under normal cir- 
cumstances, stride to stride regulation of gait is con- 
trolled by the nervous system in a largely automatic fash- 
ion [|jlO). Gait is regulated in part by "neural centers" 
within the cortex and the spinal cord As a result 

of time-varying inputs, neural activity at different times 
are thought to be dominated by different centers, lead- 
ing to complex fluctuations in the ISI "output." Further- 
more, as the central nervous system matures from infancy 
to adulthood, the interaction between neural centers be- 
comes richer 

To understand the underlying regulatory mechanisms 
of walking, deterministic and stochastic models have 
been proposed. For example, classic "central pattern gen- 
erator" models are based on oscillatory neural activity, 
where the interaction between neural centers helps reg- 
ulate gait dynamics A stochastic version of a cen- 
tral pattern generator model reproduces certain fractal 
properties of the ISI series ||^. However, existing mod- 
els do not explain observed changes in scaling exponents 
1^ , and volatility (magnitude) correlations (ijjl^ (a new 
flnding reported in the present study) that occur during 
gait maturation from childhood to adulthood. 

We propose a stochastic model consisting of a random 
walk (RW) on a chain, the elements of which represent 
excitable neural centers . A step of the RW between 
element i and element j represents the "hopping" of the 
excitation from center i to center j. The increase of neu- 



ral interconnectedness with maturation is modeled by in- 
creasing the range of "jump" sizes of the RW, since larger 
jump sizes will allow exploration of more neural centers. 
This property mimics one aspect of the increasing com- 
plexity of the adult nervous system. 

Previous studies have identified neural centers 
with pacemaker-like qualities that fire with frequency /i, 
so we represent the network of neural centers by different 
frequency modes. One mode is activated at a given time 
{ISI cx l//i), and the /, are Gaussian distributed. The 
model is based on the following assumptions (Fig. |^) : 

• Assumption (i) is that the fi have finite-size corre- 



lations, {fifi+s)/{ft) 
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We assume finite-range 



correlations among fi because neighboring neurons are 
likely to be influenced by similar factors |]lO|. This as- 
sumption effectively creates "neuronal zones" composed 
of neural centers (modes) along the chain with a typical 
size ^Q. 

• Assumption (ii) concerns the rule followed by the RW 
process. The active neural center is determined by the 
location of the RW. The "jump" sizes of the RW follow 
a Gaussian distribution of width C. 

• Assumption (iii) is that a small fraction of noise is 
added to the output of each mode to mimic biological 
noise not otherwise modeled. The output y becomes 
y{l + All) where A is the noise level and rj is Gaussian 
white noise with zero mean and unit variance [Il6|. 




FIG. 1. Illustration of the "neural" hopping model. The 
values of fi are not uncorrelated but rather have a finite size 
correlations. Shown is a sequence of four transitions, from 
mode /i+4 to fi+7 to /i+3 to fi to fi+2 ■■■ Larger values of the 
hopping-range parameter C axe associated with larger "jump 
sizes" along the chain. The neuronal zone of size 5o = 4 is 
indicated by the dashed boxes. 
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FIG. 2. (a) Examples of ISI series of healthy subjects, ages 
5 and 25 years, (b) Examples of ISI generated by the model. 
Iterating the model with a small value of the hopping-range 
parameter (C = 3) mimics the ISI of young children, while a 
large value (C — 25) mimics that of adults. 

The model has three parameters Sq, C, and A. We find 
that the best agreement virith the data is achieved when 
A = 0.02 and Sq = 25. In order to simulate changes with 
maturation, we vary only the third parameter, C, as a 
function of age, C = (age - 2) for ages 3 to 25 years. 
Increasing the hopping range with age is consistent with 
the fact that neural transmission is not fully developed 
until the late teens P,p7[|. 

Examples of ISI time series are shown in Fig. The 
ISI series of the adult subject has smaller fluctuations 
compared with the more variable ISI of the child. Two 
examples of the model's output are shown in Fig. ^d; 
using a large value of C (C = 25) simulates the ISI series 
of an adult, while using a small value (C = 3) mimics 
the ISI series of a child. Surprisingly, we find that the 
normalized probability distributions of the child ISI in- 
crement series, A{ISI)i = {ISI)i+i — {ISI)i, and of the 
model with small C are close to exponential (Fig. |3^). 
In contrast, the distribution converges to Gaussian in the 
adult and for large C in the model (Fig. ||b). 

To further test the model, we study time correlations 
in the ISI series. We calculate the function F{n), which 
corresponds to the rms fluctuations of the integrated ISI 
series of the adult (Fig. |2^) and of the model for C = 25 
(Fig. ||b) [HHjl^ . Estimating the scaling exponent from 
the slope of logF(n) vs logn, we find, for both data and 
model, that the ISI time series have correlations with 
scaling exponents ~ 0.8 (Fig. |^a). 

The scaling exponents of the original signal provide 
an indication of the linear properties (two-point correla- 
tions). Certain nonlinear aspects may be associated with 
the presence of long-range correlations in the magnitudes 
of ISI increments |A(/S'/)i|, an index of "volatility" [^ . 
We find (Fig. IJb) that the adult magnitude series is un- 
correlated (the scaling exponent is close to 0.5). The 
model shows similar behavior for large C (Fig. Qd). 



Next we inquire if changes in gait dynamics from child- 
hood to adulthood might be reflected in changes in the 
scaling exponents of the ISI series. The short-range scal- 
ing exponents of the ISI series decrease as children ma- 
ture Q. We compare the short range scaling exponents 
for a group of 50 children ^ with those of 10 adults [|; 
this exponent decreases from ~ 1.0 to ~ 0.7 (Fig. |^). By 
altering C, the model simulates these maturation-related 
changes in two-point correlation properties. 
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FIG. 3. (a) The normalized distribution of the increments 
of ISI series (A(/5'/)i = (/S'7)i+i - {ISI)i) for the child's 
ISI series (Fig. ^) and for the model with C = 3 (Fig. 
^); both model and data are consistent with an exponential 
distribution, P{x) = e~^'^^ (solid line), (b) Same as (a) for 
the adult shown in Fig. ^ and for model (C — 25) shown in 
Fig. in this histogram, both data and model are consistent 
with a Gaussian distribution, P{x) = e~'^^ (solid line). 
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FIG. 4. (a) The rms fluctuation function F{n) for the ISI 
series of adult (open circles) and of the simulation (C = 25, 
black squares) shown in Fig. |^ Here n indicates the win- 
dow size in stride number. Both the data and the model 
have similar long-range correlation properties for the ISI se- 
ries, (b) In contrast, F(n) of the magnitudes of the ISI incre- 
ments I A (75*/) i I shows weak correlations (scaling exponent of 
~ 0.55) where the scaling exponent of an uncorrelated series 
is fa 0.5. 

The magnitude series exponent of the ISI series also 
decreases with maturation (Fig. ^p). Our results for the 



2 



magnitude series suggest that the gait pattern of children 
is more volatile (and thus more nonlinear) than the 
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FIG. 5. The scaling; exponents of the ISI series of the gait 
maturation database |1£|| and for the model. The gait mat- 
uration database consists of ISI time series from 50 children 
between 3 and 13 years old. Each time series is around 8 min- 
utes (~ 500 data points). To study the effects of maturation, 
we divide the database into 5 subgroups: (i) 3-4 year olds (11 
subjects), (ii) 5-6 year olds (10 subjects), (iii) 7-8 year olds 
(14 subjects), (iv) 10-11 year olds (10 subjects), and (v) 12-13 
year olds (5 subjects). We also show data fel from an adult 
group (10 subjects 1 hour long each; ages 20-30 years). Values 
for the scaling exponent axis are given as mean ± standard 
deviation. For the model simulation, we generate 40 realiza- 
tions for each value of C; the average value is presented. The 
standard deviation around the average is ~ 0.06 for the orig- 
inal series and ~ 0.04 for the magnitude series. The age axis 
for the model follows the relation: age (years) — C + 2. (a) 
The short-range scaling exponents of the original time series 
both for the data (open circles) and the model (black squares) . 
The exponents calculated for window sizes 6 < n < 13 steps, 
decrease with age [^. The scaling exponent obtained by the 
model decreases monotonically as C increases and is within 
the error bars of the data, (b) The scaling exponent of the 
ISI magnitude series, \ A{ISI)i\. The magnitude scaling expo- 
nent is calculated for window size 6 < n < 64 for the children's 
group and 6 < n < 256 for the adult's group; the maximum 
window size is ~ .jij of the series length for both groups. The 
magnitude scaling exponent decreases with age, indicating a 
loss of magnitude correlations with maturation. The model 
exhibits a similar decrease and the simulation is within the er- 
ror bars of the actual data. The subject-to-subject variability 
is consistent with the scatter observed in physiologic indices 
of neural development . 

usual walking pattern of adults. We note that: (i) un- 
like adults, young children have difficulty keeping their 
walking speed constant. A large "walking error" in one 
direction is likely to be followed by a large compensatory 
walking error in the opposite direction — i.e., a large (or 



small) increment of ISI is likely to be followed by a large 
(or small) decrement over a given range of scales. This 
instability may lead to increased magnitude correlations, 
(ii) Adults can voluntarily simulate the less stable gait of 
children and thereby increase the magnitude correlations 
of the ISI. Children, however, cannot maintain the less 
volatile dynamics of adults. Thus, in general, the more 
mature adult gait dynamics are likely to be richer. The 
model shows (Fig. ||b) a similar decrease of the magni- 
tude exponents when increasing C and is within the error 
bars of the data. 
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FIG. 6. The multifractal formalism applied to time se- 
ries generated by the model. For the multifractal analysis, 
we applied the wavelet transform modulus maxima technique 
po[ with the Daubechies 10-tap discrete wavelet |Q. We 
integrate the series before applying the multifractal formal- 
ism. The parameters of the model are fixed (as described 
in the text) and only C changes. The series length is 32768 
data points and for each value of C, we generate 10 realiza- 
tions; the average value is presented. As the first step of the 
multifractal analysis we extract the partition function Zq{n) 
as a function of the window size n for different moments q. 
Then we estimate the scaling exponent r(g) from the slope of 
logj Zq{n) vs log2 n where 4 < n < 2048. (a) The exponent 
r(g) vs the moment q for C — 2. (b) The Legendre trans- 
form D{h) = qh — T{q) of (a) where h{q) = dr/dq and D{h) 
is the multifractal spectrum. The maximum and minimum 
Holder exponents, /imax and /imin are indicated by dashed 
lines and estimated as the slope of r(q) for —20 < g < — 15 
and 15 < g < 20, respectively (see the dashed lines in (a)). 
Error bars (1 standard deviation) are shown for both axes for 
q — —20, g = 0, and g — 20. (c) /imax and /imin for different 
values of C. A smaller C value yields a broader multifrac- 
tal spectrum. The error bars are as indicated in (b) . (d) The 
broadness of the multifractal spectrum, femax — ftmin, as a func- 
tion of C. For larger C, the multifractal spectrum becomes 
narrower. For C ^ 1 ftmax ^ /imin ~ 0.5. 

The decrease of the magnitude exponent with matura- 
tion is consistent with the possibility that the ISI time 
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series (during goal-directed walking) becomes less multi- 
fractal as the individual matures. This prediction cannot 
be tested on the available data since the time series of 
the children are too short for multifractal analysis; but 
this prediction can be tested on the model. Fig. || sum- 
marizes a systematic multifractal analysis on ISI series 
generated by the model. We find narrower multifrac- 
tal spectra with increasing C, suggesting the possibility 
of a decrease of multifractality with maturation during 
normal walking. The model's prediction agrees with a re- 
cent finding regarding the monofractality of normal gait 
in adults The model also serves as a generator of 
random series with different multifractal spectra — mul- 
tifractal for C 1 and monofractal for C ^ 1. 

In summary, we find that a simple stochastic model 
captures multiple aspects of gait dynamics, and their 
changes with maturation, including: (i) the shape of the 
probability distribution of the ISI increments; (ii) correla- 
tion properties of the ISI; and (iii) correlations properties 
in the magnitudes of the ISI increments ||2^. Further, by 
varying only a single "hopping-range" parameter, C, a 
wide array of multifractal dynamics can be generated. 
The model can also be altered by "knocking out" certain 
frequency modes (akin to what may occur during very ad- 
vanced age or in response to neurodegenerative disease). 
Simulation with drop-out of frequency modes predicts 
increased gait variability, with (i) increased magnitude 
exponents, and (ii) decrease of long-range correlations. 
Our preliminary analysis of the ISI series of older adults 
prone to falls is consistent with this prediction. Gener- 
alization of the model to two and three dimensional net- 
works to describe other types of neurological activities is 
underway. 
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